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Abstract. Pulsation powered by deuterium burning in brown dwarfs and very low mass stars has been put forth (Palla & 
Baraffe 2005) as a novel probe of the interiors of these objects in the 1-15 Myr age range. Previous observations have hinted 
at variability on the expected timescales of a few hours, suggesting but not confirming that the phenomenon is at work in 
young brown dwarfs. We have recently carried out a dedicated campaign to search for this putative class of pulsators among 
known low-mass members of five young star clusters. Our survey achieves sensitivity to periodic oscillations with photometric 
amplitudes down to several millimagnitudes. We present the census of variability over timescales ranging from minutes to 
days and discuss the current prospects for pulsation as a tool in the study of young, objects near the substellar boundary. As a 
byproduct, this work provides new insights into the distribution of stellar rotation periods at young ages via the detection of 
variability due to cool surface spots. 
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PULSATION AS A WINDOW INTO 
YOUNG BROWN DWARFS 

With masses between approximately 0.013 and 0.08 M e 
(13.6-83.8 Mj up ), brown dwarfs (BDs) provide a crucial 
link between stellar and planetary evolution. Theoretical 
models [Ql 0] are moderately successful in reproducing 
the observed properties of BDs, at least at intermediate 
and old ages (3J. Such models can be used to estimate 
fundamental parameters, namely mass and age. However 
the models have not been extensively tested by indepen- 
dent measurements of these fundamental quantities, nor 
do they account in detail for more complex realities such 
as magnetic fields and rotation which can have signif- 
icant effects as illustrated for very low mass stars by 
0]. Over the past dozen years, substantial populations 
of brown dwarfs have been identified at relatively young 
ages and in star forming regions (e.g., |g]; 0]; |§]). 
They are observed to host intriguing phenomena such as 
disk accretion and mass outflow, similar to their more 
massive brethren, in addition to magnetic activity and 
atmospheric weather, common at all ages among BDs. 
What is needed now is correlation of theoretical predic- 
tions with observational phenomena in order to better 
constrain - using fundamental physics - the basic param- 
eters of brown dwarfs. We discuss here one such possible 
new method. 

Brown dwarfs are not typical targets for stellar seis- 
mology. However, the wealth of recently discovered 
young cluster populations as well as new theoretical 
modeling efforts offer incentive to search for pulsations 
in these objects. Palla & Baraffe ([9]; hereafter PB05) 



have proposed that young (~1-15 Myr) BDs and very 
low mass stars (VLMSs) are subject to radial instability 
as a result of central deuterium burning and its strong 
temperature dependence (e.g., the e mechanism). Based 
on non-adiabatic models, they predict significant mode 
growth in a mass range from the BDs (M^0.02 M©) to 
the VLMSs (M~0.1 M ). The corresponding pulsation 
periods are expected to lie between one and five hours. 
Using a different set of models from [1] and a linear pul- 
sation code, we have independently verified these basic 
results. Even more encouraging is the substantial over- 
lap between available observational data for young BDs 
and VLMSs and the pulsation instability strip proposed 
by PB05, presented in Fig. 1. 

Although theoretical predictions provided the initial 
motivation to search for pulsation in young, low-mass 
cluster members, some observational evidence of pe- 
riodic variability has been reported on the appropriate 
^hour timescales. Several particularly promising pulsa- 
tion candidates have been uncovered in the a Orionis 
cluster llOn . II ill , and III 2|| . Such variability is unlikely 
to be the result of rotational modulation of surface spots, 
since the short periods are inconsistent with rotational 
velocities previously derived for a handful of these ob- 
jects m. 

We present here a dedicated observational campaign 
to photometrically monitor previously identified low- 
mass cluster members. Our observations are aimed at 
confirming the phenomenon of pulsation in young BDs, 
and could lead to the first anchoring of the low-mass end 
of interior models - upon which all studies of the initial 
mass function and star formation history are based. 



OBSERVATIONAL CAMPAIGN 



Search Space 

We have set out to perform a comprehensive photo- 
metric survey of BDs and VLMSs in a number of 1- 
5 Myr star clusters. At these ages, objects below 0.1 M© 
experience the deuterium burning stage of the pre-main 
sequence, making them potential pulsation candidates. 
We assembled a target sample from a set of confirmed 
and likely members of the IC 348 cluster, Upper Scor- 
pius, Taurus, Chamaeleon I, and the a Orionis cluster. 
Spectral types are available for many objects, and we 
have employed the calibration of IU4I1 to convert them 
to the temperature scale underpinning the brown dwarf 
models in PB05 [2]. With a combination of photomet- 
ric and spectroscopic data, we are able to place most 
known very low-mass young cluster members on the H- 
R diagram. The result, as shown in Fig. 1, indicates a 
promising number of BDs and VLMSs on or near the 
instability strip. Because of measurement errors and pos- 
sible model-dependent systematic effects in the position 
of the strip, we have sought to survey a large swath of 
parameter space above the zero-age main sequence. 

Not only must our survey target objects favorable for 
pulsation, but its design must also take into account ex- 
pected variability properties and limits. The linear the- 
ory of pulsation in BDs and VLMSs predicts periods 
of ~l-4 hours, but cannot determine photometric vari- 
ability amplitudes or the effect of convection on pulsa- 
tion. We therefore designed observations to probe for 
pulsation to several millimagnitudes on timescales of a 
few hours. Since sensitivity to periodic signals depends 
on both signal-to-noise (S/N) and total number of data 
points, the detection limits reached with a 1 -meter class 
telescope observing over one or two weeks are roughly 
equivalent to those achieved with a few nights with a 
larger telescope. We have selected the former set-up, 
since the longer time baseline offers increased frequency 
resolution and lower probability of weather interruptions. 

The majority of our observations were carried out with 
the CTIO 1.0-meter telescope and the robotic Palomar 
60-inch telescope. For each selected cluster field, we ob- 
serve for a total duration of 10-20 nights with continu- 
ous monitoring at 5-15 minute cadence (apart from oc- 
casional weather interruptions). We have aimed for 1% 
or better photometric precision, enabling identification 
of periodic variability in BDs with amplitudes as low 
as 0.005 magnitudes. Data are primarily in the / band, 
where late-type objects are relatively bright but light is 
not dominated by potential dust disks. 

These observations provide an unprecedented combi- 
nation of cadence, precision, and comprehensive cover- 
age of very low mass cluster members. As such, we have 
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FIGURE 1. Many young cluster members are now known to 
fall on or near the deuterium-burning pulsation instability strip 
from PB05 (hashed area). Plotted here are effective tempera- 
tures and luminosities (in units of L©) for BDs and VLMSs 
in the five clusters included in this photometric survey. As ex- 
pected from the ~l-5 Myr age estimates, all points are well 
above the zero age main sequence (ZAMS). Data were drawn 
primarily from Ql, (TJ], Q3], (Tj, (H, and d] and have 
typical errors of 100 K in temperature, and 0.15 dex in log 
luminosity. 



set out to confirm or place limits on the existence of pul- 
sation by covering a statistically significant sample of 
deuterium-burning objects across different masses and 
ages. The campaign is nearing completion, with obser- 
vations thus far of most BDs and VLMSs in IC 348, a 
Orionis, and Upper Scorpius, as well as about half of the 
intended survey regions in Cha I. Table 1 lists details on 
the five young clusters we have targeted, as well as ob- 
servational time acquired; numbers of BDs and VLMSs 
listed refer to the entire cluster, and we list the fraction 
of these monitored. 



Variability Detection 

We have produced lightcurves for data obtained thus 
far with a combination of aperture and image-subtraction 
photometry. In non-crowded fields, we employ a mod- 
ified version of the variable-aperture photometry pro- 
gram Vaphot II20I1 . Initial aperture sizes are chosen on an 
object-by-object basis so as to maximize the S/N ratio. 
Aperture sizes are then scaled with the seeing full-width- 
half-maximum (FWHM) in subsequent images. This ap- 
proach conserves the fraction of stellar flux measured 
while maintaining optimal S/N. For objects with close 
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companions whose light contaminates the aperture or sky 
background, we have instead relied on the Hotpants im- 
age subtraction package I2TI1 to remove flux from all non- 
variable objects. We then perform aperture photometry 
on the output residual images. These procedures have re- 
sulted in photometric performance very close to the lim- 
its imposed by poisson noise and sky background. The 
typical lightcurve RMS as a function of /-band magni- 
tude is presented in Fig. 2. The resulting light curves dis- 
play several different types of variable behavior which 
we have investigated for signs of pulsation. 

According to the theory of PB05, pulsation in BDs 
and VLMSs is expected to induce brightness fluctuations 
at a single frequency, which in theory can be identified 
with a power spectrum. However, it is not known a priori 
how many of the objects in our sample are host to mul- 
tiple types of variability, some of which may be high- 
amplitude and non-periodic. The plot of object magni- 
tude versus lightcurve RMS in Fig. 2 offers a quick vi- 
sual assessment of the global variability properties of our 
BDs and VLMSs. Objects whose lightcurves are incon- 
sistent with a constant source appear as outliers, and we 
can statistically select them using a % 2 test. We model the 
median measurement uncertainty as a function of magni- 
tude with the estimated poisson noise, sky background, 
and a small systematic error component. For variability 
detection at >99% certainty, lightcurves must have a % 2 
value greater than 6.6. Equivalently, we identify all data 
points lying more than a factor of 2.58 above the median 
photometric uncertainty as very likely variables. 

To probe lightcurves for lower-level variability and 
search specifically for periodic signals, we employed the 
program Period04 I22I1 . which computes an oversampled 
discrete fourier transform and may also be applied to 
time series with gaps. Periodic lightcurve variations ap- 
pear as peaks in this periodogram. Since variability is 
concentrated in a narrow range of frequency space, we 
are able to probe to amplitudes at or below the level of 
the photometric precision, particularly at periods of <5 
hours where noise is mostly uncorrelated. Based on ll23ll . 
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FIGURE 2. Example of photometric performance as a func- 
tion of magnitude, from our data on the o" Orionis cluster. 
Dots represent photometry of the entire field, while boxes de- 
note likely cluster members and crosses mark objects that have 
passed one or more of our variability detection criteria. 



we have selected a criterion of S/N>4 in the periodogram 
(where signal is the amplitude of variability) for detec- 
tion of periodic variables at 99% confidence. For typical 
BDs with / ~ 17 and an expected pulsation period of \-A 
hours, we are sensitive to variability amplitudes down to 
several millimagnitudes, as seen from the periodogram 
noise level in Fig. 4. 

RESULTS 

Rotation and Accretion-related Variability 

While this study was initially designed to search for 
the putative new effect of pulsation in deuterium-burning 
objects, the large dataset has also uncovered numerous 
other phenomena, such as brightness modulation by ro- 
tating stellar spots, accretion-related fluctuations, and 
eclipses by companions. These effects are present in ob- 



jects across a wider range in mass than is expected for 
pulsation. Thus we can analyze and classify variable be- 
havior in slightly higher mass stars (M ~ 0.1 — 0.5Mq), 
and use the results to help interpret lightcurves obtained 
for the subsample of BDs and VLMSs. 

We find that within the M0-M9 spectral type range 
(e.g., M ~ 0.02 - O.5M ), fully 50% of young low-mass 
cluster members surveyed are variable, as evident in Fig. 
2. Approximately one third of these exhibit variations 
at the — 10% level but lack detectable periodicities. As 
shown in the example at the top of Fig. 3, such aperiodic 
objects are likely accretors. In contrast, most of the re- 
maining variables are unambiguously periodic with am- 
plitudes at the few percent level and likely indicative of 
rotating spots, as seen at the bottom of Fig. 3. Their pe- 
riod distribution peaks near 1.8 days, with a steep drop- 
off at 15 hours and a tail toward the observing duration of 
10-15 days. The average amplitude of periodic variables 
is 0.03 magnitudes, although this value may be biased by 
our sensitivity limits. 



Signs of Pulsation in Brown Dwarfs? 

Periodic variability at the few percent level is cer- 
tainly common amongst low-mass pre main sequence 
cluster members. But can any of it be attributed to pulsa- 
tion instead of spots? In the absence of multi-band data, 
timescale is the main factor available to distinguish the 
two phenomena. While rotation periods of less than one 
day are not unheard of for young very low mass stars, 
variability on 1-2-hour timescales is exceedingly short to 
be explained by the passage of rotating spots. Typical ro- 
tation periods seen amongst young cluster objects seem 
to fall around 1-10 days |24J], although there is some ev- 
idence that this decreases into the brown dwarf regime 
rf25ll . In addition, rotational velocities corresponding to a 
1-hour period would reach 100 km s _1 and higher, un- 
comfortably close to the break-up speeds estimated for 
objects at these ages. Consequently, we consider any ob- 
jects with variability periods under —5 hours to be prime 
pulsation candidates. 

Despite sensitivity to millimagnitude amplitudes, the 
results of our periodogram analysis indicate a dearth of 
variables with periods less than 0.7 day. Furthermore, 
most of the brown dwarf targets show a surprising lack of 
variability on any timescale. Contrary to previous reports 
C lfToll ; fiHl ; lfl2ll ). we do not observe definitive variabil- 
ity with periods shorter than 15 hours. This is somewhat 
unexpected given the copious spot-related variability ob- 
served in slightly higher mass objects, and may be in- 
dicative of physical trends in the surface properties of 
brown dwarfs. 

Interpretation of these results depends crucially on 
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FIGURE 3. Lightcurves for two very-low-mass stars in 
IC 348 and o~ Orionis. Top: High-amplitude variability in this 
VLMS was confirmed with the % 2 test, whereas power spec- 
trum analysis did not reveal any periodic frequencies. The 
stochastic nature of the variability thus suggests that the ob- 
ject is actively accreting. Bottom: A typical periodic variable, 
with 0.03-magnitude amplitude. The lightcurve fluctuates with 
a period of 42 hours, a timescale suggestive of rotating spot 
features. 



how many objects in the sample are actually on the 
narrow instability strip. Because of field-of-view con- 
straints, we have not had the luxury of specifically tar- 
geting many of the best pulsation candidates. But obser- 
vational uncertainties (typically —0.15 dex in log lumi- 
nosity and 100 K in temperature), possible intracluster 
age spreads, and the uncertain location of the strip itself 
also create difficulty in determining which objects are the 
best pulsation targets. Assuming that the width of the in- 
stability strip is approximately correct as calculated by 
PB05, we estimate that less than 30% of the BDs and 
VLMSs in our photometric sample falls on or near it. If 
we count only objects that are directly on the strip, this 
fraction drops to 22%. Out of -50 BDs and VLMSs in 
the sample, we therefore night expect —10-15 pulsators. 

In reality, we have identified just a handful of cases 
thus far in which possible periodicities appear in the pe- 
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FIGURE 4. Top: The periodogram for a young brown dwarf 
in our sample shows a potential periodicity (marked with the 
dashed line) with amplitude 0.0044 magnitudes at 20.1 cy- 
cles/day (1.2 hours). At low frequencies the data suffer from 
correlated ("red") noise, but beyond 10 cycles day -1 the noise 
level in the periodogram reaches 0.001 magnitudes. Hence this 
tentative detection appears at signal-to-noise of just over 4.0. 
Bottom: The light curve is phased to the possible signal. Due 
to the marginal nature of the detection, the signal cannot be 
confirmed by eye. 



riodogram at 1-5-hour timescales. All of these are just 
above the detection limit at S/N;>4.0; an example is pre- 
sented in Fig. 4. Pulsation may be present in these 4-5 
objects, but it is at such a low level that confirmation with 
a different instrument set-up is required. Future observa- 
tional plans include follow-up observations of some of 
these most promising targets with both the Spitzer Space 
Telescope and Hubble Space Telescope. 

In conclusion, we have presented an extensive pho- 
tometric monitoring campaign on 1-5 Myr BDs and 
VLMSs, in hopes of confirming whether some of these 
objects are host to a new class of pulsation. While current 
evidence is not definitive, upcoming space-based data 
may provide the final word on variability in young, very- 
low-mass objects. 
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